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Phase Diagram of a 2D Vertex Model 

Hiroshi Takasaki, Tomotoshi Nishino and Yasuhiro Hieida 

Department of Physics, Faculty of Science, Kobe University, Rokkodai 657-8501 
(Received February 1, 2008) 

KEYWORDS: Vertex Model, DMRG, CTMRG, Phase Transition 



The two-dimensional (2D) vertex model is one of the 
most intensively studied systems in statistical physics. a 
Analytic expressions of thermodynamic and correlation 
functions are known in the jparameter area where Yang- 
Baxter relation is satisfiedlm- 1 The phase diagram out- 
side the solvable area is, however, not known completely. 

As a small step to explore the phase diagram of the 
16-vertex model, we investigate its subset which allows 
7 vertex configurations shown in Fig. 1, where thick and 
thin lines, respectively, correspond to the spin states 
and 1 on the bond.tr The Yang-Baxter relation is not sat- 
isfied since there is no configuration where two Is touch 
such as (1100) and (0110). 

If we interpret and 1, respectively, as absence and 
presence of adsorbed atoms on the solid surface, the con- 
dition in Fig. 1 is equivalent to setting the minimum dis- 
tance between adsorbed atoms as the lattice constant. In 
other .word, the model is equivalent to the hard square 
modelEf on the diagonal lattice. 



(0000) 



(1010) 



(0101) 



(1000) 



(0100) (0010) 



(0001) 



Fig. 1. Allowed vertex configurations, that are (abed) = (0000), 
(1010), (0101), (1000), (0100), (0010), and (0001). 



An alternative interpretation of the model is, as it has 
been done for the 6-vertex model, to regard it as 2D 
lattice polymer on the solid surface, where (1) rep- 
resents absence (presence) of a unit monomer on the 
bond. According to the allowed configurations in Fig. 1, 
only straight polymers are allowed to exist. As shown in 
Fig. 2, it is easily imagined that ordered phase appears 
when the polymer length is long or its density is high, 
and that disordered phase appears under the opposite 
conditions. The aim of this short note is to draw the 
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phase boundary of this order-disorder transition. 



disordered 



ordered 



Fie 



Disordered state and the ordered state. 



Let us parameterize the vertex model from the view 
point of the lattice polymer. For simplicity, we assume 
that the local Boltzmann weights are symmetric under 
the right angle rotations. Each vertex weight is then 
expressed as 

1^(0000) = exp(0) = 1 

w(ioio) = w(oioi) = cx P (a:) 

VF(1000) = W(0100) = VF(0010) (1) 

= W(0001) = exp(iv72 + B) 

using two parameters K = —(3U and B = —f3V, where 
U is the energy per length of the line polymer and V is 
its boundary energy; the Boltzmann weight for a straight 
line of length N is exp(NK + 2B). 

An appropriate order parameter for the order-disorder 
transition in Fig. 2 is 



M = P(1010) 
- P(0101) 



P(1000) + P(0010) 



P(0100) + P(0001) 



(2) 



where P(abcd) is the absolute probability to observe the 
local configuration (abed). We calculate M at the cen- 
ter of square systems with linear dimension L using the 
CTMRG method,© which is a variant-qf the density ma- 
trix renormlaizatipn group (DMRGp'tP applied to 2D 
classical systemSjEf where L is chosen to be sufficiently 
larger than the correlation length. To weakly stabilize 
the order in the horizontal direction, we impose fixed 
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boundary conditions where all the spins on the vertical 
sides of the square are 1 and those on the horizontal sides 
are 0. This boundary condition makes the corner trans- 
fer matrix asymmetric, therefore we use an asymmetric 
extension of the CTMRG method. 



003 



1.334 1.336 1.338 1.340 1.342 1.344 1.346 

K 

Fig. 3. The eighth power of the order parameter (= M 8 ) with 
respect to K when B = 0. 



The phase boundary is determined from the observa- 
tion of the order parameter M . Figure 3 shows the in- 
dependence of M 8 when B = 0. It is clear that M 8 is 
a linear function of K, and we obtain the critical point 
K c = 1.334. Also for the cases where B / we have 
checked that M 8 is linear in K; the critical index (3 of 
this system is 1/8. 




Fig. 4. Phase diagram in the K-B plane. 



We thus obtained the phase boundary in the K-B 
plane shown in Fig. 4. The system is always disordered 
in the negative K region, therefore we only show the 
positive K region. The boundary starts from (K, B) = 
(0 + ,-2.25), it passes (K, B) = (1.334,0), and then ap- 
proaches the line B = (K — 1.0)/2; we have not iden- 
tified the reason why the line has the offset of amount 
-1.0/2 = -0.5. 

It is worth obtaining another critical index of this 



model, which is related to the density defined by 

P(1000) + P(0010) 



n = P(1010) + 
+ P(0101) + 



P(0100) + P(0001) 



(3) 



Figure 5 shows the K dependence of n — n c when B = 0, 
where n c = 0.736 is the density at the phase boundary 
K c = 1.334. It is apparent that n — n c is proportional to 
K — K c . Since n is related to the energy parameter U, 
the exponent a is equal to 0. 



n-ric 



0.000 ' 1 1 1 

0.000 0.005 0.010 0.015 0.020 



K-Kc 



Fig. 5. Density in Eq. 3 near the transition point K c = 1.334 
when B = 0, where n r = 0.736. 



In conclusion we have studied a symmetric vertex 
model which allows 7 vertex configurations shown in 
Fig. 1, and obtained its phase diagram. The critical in- 
dices of this model, that are (3 = 1/8 and a = 0, are 
the same as those of the Ising model. To include ad- 
ditional configurations, such as (1111) and (1110), that 
may cause the percolation transition, is our next sub- 
ject toward the complete understanding of the 16-vertex 
model. 
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